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Abstract - In pursuit of rational control of orbital polarization, we present a combined experi- 
mental and theoretical study of single unit cell superlattices of the correlated metal LaNiC>3 and 
the band insulator LaAlOs. Polarized x-ray absorption spectra show a distinct asymmetry in the 
orbital response under strain. A splitting of orbital energies consistent with octahedral distortions 
is found for the case of compressive strain. In sharp contrast, for tensile strain, no splitting is 
found although a strong orbital polarization is present. Density functional theory calculations 
including a Hubbard U term reveal that this asymmetry is a result of the interplay of strain and 
confinement induces octahedral rotations and distortions and altered covalency in the bonding 
across the interfacial Ni-O-Al apical oxygen, leading to a charge disporportionation at the Ni sites 
for tensile strain. 
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The interplay between structure and strong correlations 
in transition metal (TM) oxides gives rise to a breadth of 
intriguing phenomena ranging from high-temperature su- 
perconductivity to magnetism [1] . The ground state prop- 
erties are fundamentally influenced by the orbital occupa- 
tion. For the case of charge and spin degrees of freedom, 
electric and magnetic fields provide a direct means to con- 
trol the energy of those elements of the ground state. For 
the case of orbital degrees of freedom the route to direct 
control is not well established. With the recent advances in 
oxide film synthesis, strain fields in pcrovskitc hctcrostruc- 
turcs are by far the most common route used to attempt 
rational control of orbital energies [Tj. However, recent 
experiments have demonstrated that altered chemical co- 
ordination at the interface can lead to a reconstruction of 
the orbital configuration into a unique correlated electron 
states not found in the bulk [JHSJ. 

It was recently theoretically proposed that through 
a combination of strain and spatial confinement one 
can possibly manipulate the orbital state in heterostruc- 
tured LaNiOa (LNO). In bulk, LaM0 3 , the Ni 3+ (3d 7 , 
t2 g e*,S = 1/2) forms a paramagnetic metallic state [7] 
with no evidence of orbital preference |S|, but it was the- 
oretically suggested that strain and confinement could 



break this symmetry in LaNi0 3 /LaM0 3 (M=Al,Ga,...) 
superlattices and thereby allow the manipulation of the 
orbital occupancy [S]. Calculations within the local den- 
sity approximation in combination with dynamical mean 
theory supported the formation of a cuprate-like Fermi 
surface in the LaNiC^/LaMO^ system [TU]. Additional 
theoretical work also demonstrated that the chemical con- 
finement is more important than strain in controlling the 
orbital polarization [llj . On the experimental side, several 
recent studies have looked at several unit cell superlattices 
and shown a metal to insulator transition due to confine- 
ment effects [ElE], which was attributed to onset of a 
charge-ordered antiferromagnetic state at low temperature 
[14] . Polarization dependent resonant reflectivity was also 
utilized to demonstrate an enhanced orbital polarization 
at the interface with respect to the bulk of the LNO layer 
[15] . However, these studies are all limited to tensile strain 
and a full understanding requires a microscopic picture in 
the ultrathin limit that can deconvolve effects of strain vs. 
confinement. 

In this letter, we present a combined experimental and 
theoretical study on single-unit-cell LaNi0 3 /LaA103 su- 
perlattices to elucidate the possibility of cooperative chem- 
ical and strain control of orbitals. The results show an mi- 
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expected asymmetric orbital response to compressive vs. 
tensile strain. Compressive strain generates a small or- 
bital splitting consistent with expectations of elastic the- 
ory while tensile strain shows no evidence of crystal field 
splitting yet displays a signifcant d x 2_ y i orbital polariza- 
tion. Density functional theory (DFT) calculations with 
an on-site Coulomb repulsion term confirm the distinct 
electronic behavior and attribute it to the concerted ef- 
fects of broken symmetry at the interface, chemical con- 
finement, and a strain induced lattice distortion. 

Fully epitaxial [LaNi03(l u.c.)/LaA103 (1 u.c.)]20 su- 
perlattices (SL) were grown on (001) Ti02-terminated 
SrTi0 3 (STO) and LaA10 3 (LAO) single crystal sub- 
strates by pulsed laser deposition with in situ monitoring 
by RHEED using the interrupted growth method [TMIB] . 
The high-quality of the supcrlatticc was confirmed by 
AFM and high-resolution x-ray diffraction at beamline 6- 
ID-B of the Advanced Photon Source (APS) at Argonne 
National Laboratory. To determine the electronic prop- 
erties, the samples were studied with polarized x-ray ab- 
sorption in the soft x-ray regime at beamline 4-ID-C of 
the APS. The spectra were measured in bulk sensitive flu- 
orescence yield and aligned using a NiO (Ni 2+ ) standard 
measured simultaneously in the diagnostic section of the 
beamline. 

X-ray diffraction was used to gain insight into the su- 
perlattice crystal structure. In a pseudo-cubic form, the 
i?3c bulk LNO has a lattice parameter of ^3.83 A. Off- 
specular mapping of the (222) reflection (not shown) de- 
termined that the films are fully lattice matched to the 
substrate for both cases. The average c-axis lattice pa- 
rameter is 3.932 A and 3.823 A for samples on LAO and 
STO, respectively. The results for LAO are consistent 
with an out of plane expansion expected for a volume 
conserving scenario. Surprisingly, despite the large ten- 
sile strain on the STO-substrate, the SL shows almost the 
bulk LNO lattice parameter. This indicates that the su- 
perlattice structure is responding differently under tensile 
vs. compressive strain, which was seen in recent studies 
of LaNi03 films [19l[20]. The superlattice peaks are ob- 
servable with a half-order position, consistent with 2 u.c. 
periodicity of the superlattice [18]. 

To provide a more detailed insight into differences in 
the electronic properties as a function of strain we per- 
formed polarization dependent XAS at the Ni L edge . 
XAS probes the density of unoccupied states giving a pic- 
ture of not only the chemical state, but also the orbital 
occupancy. An important point to note is that due to the 
strong core hole interaction, the Ni L-edge probes the local 
Ni environment. From the polarization averaged XAS, we 
determine that the valence of Ni is bulk-like 3+ testifying 
to the absence of oxygen vacancies and lack of a signifi- 
cant amount of charge transfer at the interface (see Fig. 
Q3a)). The data are shown in comparison to bulk NdNi03 
above and below the metal-insulator transition. A strong 
multiplet feature develops on the leading edge of both the 
L3 and L2 edges, which is identified as a signature of local- 
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Fig. 1: (a) Polarization averaged Ni L edge in comparison with 
NdNiC>3 above and below the metal to insulator transition, (b) 
Polarization dependent x-ray absorption at the Ni L2 edge for 
both (b) compressive (LAO substrate) and (c) tensile (STO 
substrate) strain. The data show the expected response for 
compressive strain, but the the lack of energy splitting for ten- 
sile is quite unexpected. The differences in the orbital response 
are highlight in both cases by showing the x-ray linear dichro- 
ism (XLD). 

ized carriers in insulating charge-ordered (CO) nickelatcs 
[2"TH2"3"] . The spectroscopy illustrates directly that, at a lo- 
cal level, the carriers are localized in agreement with the 
macroscopic scale results of the transport measurements 
(not shown). This finding of a charge ordered state is also 
consistent with recent optical studies of ultrathin super- 
lattices [T4] , 

The effect of confinement and strain on orbital occu- 
pancies can be probed directly by light with tunable lin- 
ear polarization. Due to the strong overlap between the 
La M4 and Ni L3 edges, however, for the polarization de- 
pendence we focus on the behavior of the L 2 edge only. 
First, consider the case of compressive strain (LAO) as 
shown in Fig. [IJb). In-plane (ab) octahedral compression 
was shown with XRD to generate an out-of-plane (c) ex- 
pansion of the NiOg ocathedra, which leads to a lowering 
of the c-axis orbital (3z 2 — r 2 ) energy As anticipated, 
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Fig. 2: Total density of states for (LAO)i/(LNO)i for the 
case of tensile (STO) (solid black line) and compressive strain 
(LAO) (dashed blue line). 



this effect is clearly observed in the polarized XAS, which 
shows a shift of the c-axis polarized spectra to lower en- 
ergy as expected when lowering the energy of this state. 
This shift of ~100 meV is related to the energy scale for 
the crystal field splitting due this type of structural distor- 
tion [20], but it is important to note that the crystal field 
splitting is small compared to the bulk LNO bandwidth of 
several eV [7|. The polarization dependence is best seen 
in the difference of the ab- and c-axis spectra known as 
x-ray linear dichroism (XLD). Under compressive strain, 
the XLD shows a derivative like lincshape (see Fig. QJb)) 
consistent with the effect being largely an energy shift be- 
tween the ab vs. c spectra. 

For the case of tensile strain, under the assumption of a 
symmetric strain response, one would expect the opposite 
effect, but the data in Fig. [ljc) shows practically no en- 
ergy splitting between the two orbitals. Instead, there is 
a significant difference in peak intensity with polarization. 
The lack of energy shift strongly implies there is no split- 
ting of the e g band as expected for undistorted octahedra. 
On the other hand, even though there is no energy split- 
ting of the orbitals under tensile strain, the difference in 
intensity signals a significant departure from the uniform 
occupation of G?3 Z 2_ r 2 vs d x 2_ y 2 orbitals in bulk LNO. In 
this case, the XLD shows almost no intensity in the region 
of the rising edge of the peak due to the absence of an en- 
ergy shift, but has a large signal near the peak maximum 
due to the difference in e g orbital occupancies. 

To quantify the differences in strain dependent orbital 
occupancy, we integrate the area under the peak to count 
the number of unoccupied states (holes) which scales with 
the orbital occupancies. Proper counting of holes requires 
measurement of the entire L edge, which is not possible 
due to the La M-edge intrusion into the L 3 edge. How- 
ever, under the assumption that spin-orbit coupling in the 
3c? states is constant with strain, by integrating the L2 
edges we can extract the relative orbital occupancies [24] . 
For both strain states the in-plane hole counts are very 
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Fig. 3: (left) Projected (spin-averaged) density of states 
(PDOS) of Ni 3d and (right) electron density distribution of the 
unoccupied states ibetween Ep and 2.5 eV in (LAO)i/(LNO)i 
for (a) tensile (STO) and (b) compressive (LAO) strain. The 
right panels clearly show the reduction in hole density on the 
apical oxygen between Ni and Al, as well as the charge dis- 
porportionation at the Ni sites for tensile strain. 

similar indicating that the d x 2_ y 2 orbitals have a similar 
occupancy regardless of strain. For tensile strain though, 
there is a significant increase in the number of holes in the 
out-of-plance d 3z 2_ r 2 orbital states (decrease in electron 
occupancy of the d^ z 2_ r 2 band). By defining orbital polar- 
ization through the relative areas, then there is about 5% 
reduction in d 3z 2_ r 2 occupancy for tensile strain, while for 
compressive strain the orbital polarization is close to zero. 
The polarization for tensile strain agrees well with recent 
measurements using polarized x-ray reflectivity |15j . 

These experimental observations raise several important 
questions: How does one find an overall orbital polar- 
ization in the absence of crystal field splitting and why 
is there no crystal field splitting of the orbitals even un- 
der high tensile strain? To shed insight into these ques- 
tions, we performed density functional theory calcula- 
tions of (LAO)i/(LNO)i superlattices with the full po- 
tential linearized augmented plane waves method as im- 
plemented in the WIEN2k code [25]. Electronic correla- 
tions beyond the generalized gradient appoximation [21] 
of the exchange correlation potential are included within 
the LDA/GGA+C/ approach [27] using U = 4 eV and 
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J = 1 eV. For these calculations, the superlattices were 
strained laterally to the in-plane lattice parameters of STO 
and LAO and the out of plane parameter was set accord- 
ing to the experimentally obtained out-of-plane lattice pa- 
rameter. While previous DFT studies PHIE] have used 
a tetragonal setup we explored the influence of strain on 
the octahedral rotations by optimizing the internal param- 
eters in a monoclinic unit cell containing 20 atoms. Fig. [2] 
shows a finite density of states at Ef for the system un- 
der compressive strain (dashed line) . The two-dimensional 
character of these conducting sheets separated by insulat- 
ing LAO layers leads to an enhanced resistivity compared 
to bulk LNO. In contrast, for tensile strain a gap in the 
total density of states (Fig. [2]) of 0.3 eV occurs just above 
at the Fermi level (solid line, grey area) . This distinct be- 
havior around the Fermi level for tensile and compressive 
strain is consistent with the variation in resistivity of the 
two samples. 

As expected, all t^g orbitals are occupied at the Ni-sites 
and there is a significant occupation of both e g orbitals 
(see Fig. [3]). This is consistent with the large covalency 
leading to a 3d s L ground state, where L denotes a hole on 
the oxygen site, and thus the Ni valence state is closer to 
Ni 2+ rather than the expected Ni 3+ . The structural and 
electronic properties show a strongly asymmetric response 
to strain as seen in the experiment. The most striking ef- 
fect for tensile strain is a charge disproportionation on the 
Ni sublattice, Ni 3+ ->-Ni 3+,5 +Ni 3 -' 5 with a total difference 
in 3d occupation at the Nil (Ni 3-5 ) and Ni2 (Ni 3+<5 ) sites 
of ~ 0.2e. A similar tendency was recently reported for 
single-layer LNO films under tensile biaxial strain [20] , but 
in the LAO/LNO supcrlattice as a result of confinement 
the effect is strongly amplified leading to the opening of a 
gap of 0.3 eV just above the Fermi level. The charge or- 
dering goes hand in hand with bond-length disproportion- 
ation driven by a breathing mode with two nearly regular 
octahedra with in-/out-of-plane cation-anion distances of 
2.03/1.99 A and 1.91/1.92 A at the Nil and Ni2 sites, re- 
spectively [55]. This is very similar to the ground-sate of 
charge ordered NdNiOa [29] and occurs due to a reduction 
in bandwidth due to bond angle distortions [7.. In con- 
trast, for compressive strain the in-plane Ni-0 bond length 
(1.92 A) is significantly shorter than the out-of-plane dis- 
tance (2.00 A) consistent with a volume conserving tetrag- 
onal distortion. 

The influence of strain and the interface coordination in 
the superattices is seen clearly in the spatial distribution 
of holes shown in Fig. [3] (right panels). The key feature for 
both tensile and compressive strain is the reduced number 
of holes in the p z -orbital of the apical oxygen while the in- 
planc oxygen orbital occupancies are very similar to bulk 
LNO. For both strain cases there is ^6% reduction of the 
apical oxygen hole density which results from the stronger 
ionicity of the Al 3+ -0- vs. Ni 3+ -0-bond. To first order 
this effect is larger than that of strain, which gives rise 
to a chemical means of controlling the orbital occupancies 
in addition to strain as also noted by a recent theoretical 



study [TT|. We note that for tensile strain the hole at 
the apical oxygen of Nil (Ni 3 h<5 ) is almost completely 
suppressed due to the increased occupancy of the d 3z 2_ r 2 
orbital closing the path to form the 3c? 8 L state. 

Using the insight from theory, we can see clearly that 
the ground state of this system is a result of the inter- 
play of three effects: strain, confinement, and intcrfacial 
covalency. The asymmetric response to strain is reflected 
in the stabilization of a charge ordered state under ten- 
sile strain, which changes the unit cell volume and leads 
to isotropic crystal fields in the two distinct sites. On 
the other hand, under compressive strain the expansion 
of the out-of-plane Ni-0 bond length results in a crys- 
tal field splitting between the d x 2_ y 2 and d^ z 2_ r 2 orbitals 
and an energy lowering of the latter (see Fig. (3h>). This 
is analogous to previous results for compressive strain in 
thick LNO films that show quantitative agreement in the 
altered crystal field to the splitting of the orbitally re- 
solved bands at the T point [19] . By far though, the most 
profound effect is the reduction of the apical oxygen hole 
density. This leads to a reduction in the d^, z 2_ r 2 orbital 
occupancy under tensile strain even in the absence of crys- 
tal field splitting of the orbitals. For compressive strain, 
the strain induced crystal field should enhance the occu- 
pancy of the d3 Z 2_ r 2 orbital, which is not observed. This 
is due to the fact that the enhanced occupancy from strain 
is offset by the reduction due to the reduced covalency at 
the Ni-O-Al bond at the interface. 

In summary, the key aspect of this study is that con- 
trol of orbital occupancies results from a competition of 
several effects: correlation, lattice distortions, chemistry 
and strain. These results highlight that rational manipu- 
lation of orbitals requires an understanding that extends 
beyond the conventional picture of how strain alters elec- 
tronic properties. Chemical control offers new opportu- 
nities to control orbital configurations in ultra-thin oxide 
heterostructures even in the absence of strain. 
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